Pseudoscalar meson form factors and decays 

A.E. Dorokhov 
Abstract 

In this communication we discuss few topics related with modern experimental data on the 
physics of light pseudoscalar mesons. It includes the contribution of the pseudoscalar mesons to 
the muon anomalous magnetic moment (AMM), g — 2, the rare decays of light pseudoscalar mesons 
to lepton pair, the transition form factors of pseudoscalar mesons at large momentum transfer, the 
pion transversity form factor. 

Measuring the muon anomalous magnetic moment g — 2 and the rare decays of light pseudoscalar 
mesons into lepton pair P — > serve as important test of the standard model. To reduce the 
theoretical uncertainty in the standard model predictions the data on the transition form factors of 
light pseudoscalar mesons play significant role. Recently new data on behavior of these form factors 
at large momentum transfer was supplied by the BABAR collaboration. Within the nonlocal chiral 
quark model it shown how to describe these data and how the meson distribution amplitude evolves 
as a function of the dynamical quark masses and meson masses. 



I. INTRODUCTION 



Among hadrons the pion plays very special role because it has the mass much lighter 
than other hadronic states. In theory it appears as a result of spontaneous breaking of 
chiral symmetry and gains its small mass due to small masses of nonstrange current quarks. 
It seems that confinement forces do not play any significant role in formation of the pion as 
a bound state. 

Studies of the pion have long history and at the moment the processes with pion par- 
ticipation is one of standard hadronic background. From other side new highly-precise 
experiments related to rare decays or form factors at large momentum transfer serve as an 
important independent way to search for physics beyond the standard model and as test 
of the QCD dynamics. At the moment, there are some problems with the matching of the 
experimental data with the predictions of the SM. The most famous one is the discrepancy 
by three standard deviations between experiment [1] and SM theory j2j |3] for the muon 
anomalous magnetic moment (AMM). Another example is similarly large deviation between 
the recent precise experimental result on the rare 7r° decay into e~e + pair [1] and the SM 
prediction [5J. 

II. MUON G-2 

The theoretical studies of the muon AMM g — 2 (see for review [6T-HU]). the rare de- 
cays of light pseudoscalar mesons into lepton pairs [5j ITTHH] and the comparison with the 
experimental results, offer an important low-energy tests of the SM. 

The discrepancy between the present SM prediction of the muon AMM and its experi- 
mental determination pQ is (28.7 ±8.0) • 1(T 10 (3.6a) [2J. The theoretical error is dominated 
by hadronic corrections. They are the vacuum polarization (HVP) and light-by-light (LbL) 
contributions. The latter can be estimated only in model dependent way. The latest results 
are the calculations of the pseudoscalar hadronic channel contribution within the instanton 
quark model [15] and the calculations within the so called [16] Dyson-Schwinger equation 
approach. In [15] in the full kinematic dependence of the meson-two-photon vertices from 
the virtualities of the mesons and photons is taken into account. It is demonstrated that 
the effect of the full kinematic dependence in the meson-photon vertices is to reduce the 
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contribution of pseudoscalar exchanges comparing with the most of previous estimates and 
the result is a P s ' LbL = (5.85 ± 0.87) • 1(T 10 . 

The knowledge of the eta(') meson couplings to virtual photons is important for the 
calculation of the anomalous magnetic moment of the muon, being pseudoscalar exchange 
the major contribution to the hadronic light-by-light scattering. For illustration in Fig. [T] 
we present on the same plot the vertex Fp. 77 (p 2 ; 0, 0) in the timelike region p 2 < and the 
vertex Fp. 7 * 7 (p 2 ;p 2 , 0) in the spacelike region p 2 > as they look in the the instanton quark 
model (NxQM) and vector meson dominance model (VMD). These two special kinematics 
match at zero virtuality p 2 = 0. The remarkable feature of this construction is that the first 
kinematics is connected with the decay of pseudoscalar mesons into two photons at physical 
points i 7 p 77 (— M^; 0, 0) = gp 77 , while the second kinematics is relevant for the light-by-light 
contribution to the muon AMM. Thus, the part of Fig. [T]at p 2 < describes the transition of 
the pion-two-photon vertex from the physical points of meson masses to the point with zero 
virtuality, which is the edge point of the interval where the integrand is defined. In VMD 
type of models there is no such dependence on the meson virtuality. Thus, the value of this 
vertex at zero meson virtuality is the same as the value of the vertex at the physical points 
of meson masses, Fp^(p 2 = —M 2 ,;0,0) = Fp^P(0; 0, 0). However, the r\ and rf mesons 
are much heavier than the pion and such extrapolation is too crude. One can see that for rj 
and particularly for rf the difference between the values of the vertex at physical and zero 
virtuality points is large, Fp 77 (p 2 = — M 2 V ,;0,0) ^> Fp 77 (0; 0, 0). Thus, the contributions 
of the rj and rj mesons to the muon AMM evaluated in N%QM are strongly suppressed as 
compared with the VMD results that can only be considered as upper estimates of these 
contributions. 

In [16] the hadronic LbL contribution to the muon AMM using the framework of Dyson- 
Schwinger equation (DSE) with the result a^ s,qLoop = (13.6 ± 5.9) ■ 10~ 10 that is bigger 
than most previous estimates. It is not easy to understand this result, in particular, in 
view of discussion in [T7]. Moreover, it is instructive to compare the results on hadronic 
vacuum polarization contribution in the instanton [TB] and DSE [T9] based models. In [TS] 
it was demonstrated that the dominant HVP contribution to muon AMM comes from the 
dynamical quark loop and contribution of the p meson is highly suppressed, while in [19] 
the opposite conclusion was made. It is still open question what model is more adequate in 
these calculations. 
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FIG. 1. Plots of the ir°, r\ and rj vertices Fp* 77 (p 2 ; 0, 0) in the timelike region and Fp*- y * 1 (p 2 ;p 2 , 0) 
in the spacelike region in N^QM model (thick lines) and VMD model (thin lines). The points with 
error bars correspond to the physical points of the meson decays into two photons. The VMD 
curves for ir° and r\ are almost indistinguishable. 

III. RARE PION DECAY vr° e+e" 

The situation with the rare decays of the light pseudoscalar mesons into lepton pairs 
became intriguing after the KTeV E799-II experiment at FermiLab [4j in which the pion 
decay into an electron-positron pair was measured with high accuracy (R{P — > = 

T{p^in-)/v tot ) 

The standard model prediction gives [5j [13] 

R Thcor (it e+e") = (6.2 ± 0.1) • 10~ 8 , (2) 

which is 3.1a below the KTeV result 0. Other pseudoscalar mesons rare decays into lepton 
pair are also well suited as a test of the standard model (see Table). From experimental point 
of view the most interesting are the decays of 77 and 77' to muon pair. The branchings for 
these processes theoretically bounded as from below as well from the above [13] . Further 
independent experiments for 7r° — > e + e~ at WASAatCOSY [20] and for rj(rj' + l~ KLOE 
[2T] and BES III [22] and other facilities will be crucial for resolution of the problem with 
the rare leptonic decays of light pseudoscalar mesons. More precise data on the transition 
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form factors in wider region of momentum transfer are expected soon from the BABAR, 
BELLE (at large momentum transfer) and KEDR, KLOE [23] (at small momentum transfer) 
collaborations. These data would allow to make more accurate theoretical predictions. 



TABLE I. Values of the branchings R {P -> = T(P -»■ /T tot obtained in approach 0H3] 

and compared with the available experimental results. 



R 


Unitary 
bound 


CLEO+BABAR 
bound 


CLEO±BABAR 
±OBE 


With mass 
corrections 


Experiment 


R(ir° -> e+e") x 10 8 


> 4.69 


> 5.85 ±0.03 


6.23 ±0.12 


6.26 


7.49 ±0.38 @] 


R(r)^> x 10 6 


> 4.36 


< 6.60 ±0.12 


5.35 ±0.27 


4.76 


5.8 ±0.8 [211 [25] 


R(r) -»• e+e - ) x 10 9 


> 1.78 


> 4.27 ±0.02 


4.53 ±0.09 


5.19 


< 2.7- 10 4 [26] 


i?(r/ -»• x 10 7 


> 1.35 


< 1.44 ±0.01 


1.364 ±0.010 


1.24 




R(rj' -> e + e~) x 10 10 


> 0.36 


> 1.121 ±0.004 


1.182 ±0.014 


1.83 





There are quite few attempts in the literature, to explain the excess of the experimental 
data on the ir° — > e + e~ decay over the SM prediction, as a manifestation of physics beyond 
the SM. In Ref. [27], it was shown that this excess could be explained within the currently 
popular model of light dark matter involving a low mass (~ 10 MeV) vector bosons U^, which 
presumably couple to the axial-vector currents of quarks and leptons. Another possibility 
was proposed in Ref. [281 12S] , interpreting the same experimental effect as the contribution 
of the light CB-odd Higgs boson appearing in the next-to-minimal supersymmetric SM. The 
latter version perhaps excluded by modern experiments. 

IV. PSEUDOSCALAR MESON TRANSITION FORM FACTORS 

The main limitation on realistic predictions for above processes originates from the large 
distance contributions of the strong sector of the SM, where perturbative QCD does not 
work. In order to diminish the theoretical uncertainties, the use of the experimental data 
on the pion charge and transition form factors are of crucial importance. The first one, 
measured in e+e" 7r + 7r-(7) by CMD-2 [30], SND [31], KLOE [32], and BABAR [33] 
provides an estimate for the hadron vacuum polarization contribution to muon g — 2, with 
accuracy better than 1%. The second one, measured in e + e~ — > e + e~P for spacelike photons 
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by CELLO [34J, CLEO [35], and BABAR [36H38] collaborations and in e+e" ->■ P7 for 
timelike photons by the BABAR [39J collaboration, is essential to reduce the theoretical 
uncertainties in the estimates of the contributions of the hadronic light-by-light process to 
the muon g — 2 and in the estimates of the decay widths of P — > . The BABAR data 
|36[ [39] on the large momentum behavior of the form factors cause the following problems 
for their theoretical interpretation [lOj SJJ: 1) An unexpectedly slow decrease of the pion 
transition form factor at high momenta [36], 2) the qualitative difference in the behavior of 
the pion and 77,7/ form factors at high momenta [37], 3) inconsistency of the measured ratio 
of the 7],rj' form factors with the predicted one [59] . 

In Figs. [2J [3j [4] and [5] the data for the 7r°, 77, rf and r\ c transition form factors from the 
CELLO, CLEO, and BABAR collaborations are presented. In Figs. ||and|4) the CLEO 02] 
and BABAR [32] points, measured in the timelike region, are drawn at Q 2 = 14.2 GeV 2 and 
Q 2 = 112 GeV 2 , correspondingly, assuming that the spacelike and timelike asymptotics of 
the form factor are equal. It is seen from the Figs. [3] and |4j that the spacelike and timelike 
points are well conjugated. 




Q 2 (GeV 2 ) 

FIG. 2. The transition form factor 7*7 -> tt . The data are from the CELLO [34], CLEO [35] and BABAR 
[31)] Collaborations. The solid line is the nonlocal chiral quark model calculations and the dash-dot line is 
the parametrization (6). The dotted line is massless QCD asymptotic limit. 
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FIG. 3. The transition form factor 7*7 -> 77. The data are from the CELLO [M], CLEO [35] and BABAR 
[37] Collaborations. The CLEO results obtained in different r\ decay modes are averaged. The CLEO and 
BABAR points, measured in the timelike region 7* — > 777 [12] and [33], are drawn at Q 2 = 14.2 GeV 2 
and Q 2 = 112 GeV 2 , correspondingly, assuming that the spacelike and timelike asymptotics of the form 
factor are similar. The solid line is the nonlocal chiral quark model calculations and the dash-dot line is the 
parametrization (6). 

where a is the QED coupling constant, Mp is the resonance mass and F(P — > 77) is the 
two-photon partial width of the meson P. The axial anomaly predicts 
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Fprf (Q = 0, 0) 



(4) 



where fp is the meson decay constant. Under assumption of factorization, perturbative QCD 
predicts the asymptotic behavior of the Fp *(Q 2 , 0) transition form factors as Q 2 — > 00 [H] 

2/p 



F Pyr (Q 2 oo,0) 



Q 2 



(5) 



The perturbative QCD corrections to this expression at large momentum transfer are quite 
small [^5tT4"8] . The CLEO (and CELLO) collaboration parameterized their data by a formula 
similar to that proposed by Brodsky and Lepage in [H] , but with the pole mass being a free 
fitting parameter 



i^T(Q 2 ,o) 



1 



47r 2 /pl + Q 2 /Ap' 



(6) 



where the values of f P are estimated from @ and Q [35]: f v = 92.3 MeV, f v = 97.5 MeV, 
/,/ = 74.4 MeV. In Figs. [2j |3j [i] and [5] the parametrizations ^ with parameters A n = 776 
MeV, A,, = 800 MeV, = 859 MeV and A vo = 2.92 GeV are shown by dot-dashed lines. 
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FIG. 4. The transition form factor 7*7 ->■ rf. The data are from the CELLO [31], CLEO [35] and BABAR 
[37] Collaborations. The CLEO results obtained in different rf decay modes arc averaged. The CLEO and 
BABAR points, measured in the timelike region 7* — > 777 [12] and [33], are drawn at Q 2 = 14.2 GeV 2 
and Q 2 = 112 GeV 2 , correspondingly, assuming that the spacelike and timelike asymptotics of the form 
factor are similar. The solid line is the nonlocal chiral quark model calculations and the dash-dot line is the 
parametrization (6). 

We see that the parametrization (|6| describes the rf and t]q mesons form factors (Figs. 
|3J p| well at all measured momentum transfer. However, it starts to deviate from data 
for the pion and probably for the eta mesons form factors at momentum transfers squared 
larger than 10 GeV 2 (Figs. [2J [3]). Namely, the high momentum transfer data obtained by 
BABAR collaboration show evident growth of the form factor multiplied by Q 2 at large Q 2 
for the pion and probably for the eta. This is in contradiction with asymptotic formula ^ 
and unexpected from the QCD factorization approach [38£J The r/ c form factor is in good 
agreement with predictions |4"5] . 

In [501 EE] the asymptotic of the pseudoscalar meson transition form factor generalizing 
the result (|5l) was derived 



where D (k 2 ) is the nonperturbative inverse quark propagator with property D (k 2 — > 00) — > 
k 2 and <p n (x) is the meson distribution amplitude (DA) given in the nonlocal model by (see 

1 Small radiative corrections can be attributed to slight changes in the parameters Ap in (|6|). 
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FIG. 5. The transition form factor 7*7 — > r\c- The data are from the BABAR [38 Collaboration. The 
dashed line is the nonlocal model calculation and the dot-dashed line is the fit. (The normalization is 
arbitrary) 

for notations and the a representation [50j I5T] ) 



(f P (x) 



47T 2 /] 2 



PS,7T JO 



00 AT 

alj „xxLM 2 



p (xG mfi (xL, xL) + xG 0:m (xL, xL)) , 



(8) 



with x — (1 — x) and 



dxipp (x) 



The properties of the form factor at large Q 2 strongly depend on the end point behavior 
of the meson DA. In [52H51] it was assumed that the pion is almost pointlike and that its 
DA can be almost constant (flat) in its shape. In this case the integrand in (J7|) is poorly 
convergent and the 1/Q 2 behavior of the form factor is enhanced by logarithmic factor 
ln(Q 2 /M 2 ) asymptotically or in rather wide region of large Q 2 , where M q is dynamical 
(constituent) quark mass. At the same time the DA with suppressed end point behavior 
lead to the standard result 

The solid line in Fig. [2] is the pion transition form factor calculated in the nonlocal chiral 
quark model [50l I5T] with the mass parameters Mq = 135 MeV and strange quark mass 
M s = 250 well describes the BABAR data. There are other attempts to explain the BABAR 
data on the pion form factor within different factorization schemes [55H61] . 



By using singlet-octet mixing scheme one has for rj and rf meson form factors [T5 

1 



F^ 7 (m 2 ;Q 2 ,0) 



3V3 



[(5F U (m 2 v , Q 2 , 0) - 2F S (m 2 ; Q 2 , 0)) cos 0(m 2 )- 



-V2 (5F U (m 2 v , Q 2 , 0) + F s (m 2 ; Q 2 , 0)) sin I 



F, V7 (m 2 ,;Q 2 ,0) = ^= [(5F U (m 2 ,; Q 2 , 0) - 2F S (m 2 ,;Q 2 ,0)) sin#(m 2 ,) + 
W2(5F M (m 2 ,;Q 2 ,0) +F S (m 2 ,; Q 2 , 0)) cos #(m 2 ,; 



(9) 



From these expressions it is clear that in general it is not possible to separate the nonstrange 
and strange components of the form factors as it was attempted to do in [37] and then used 
in some theoretical works. 

The solid line in Figs. [3j |4j is the eta and eta prime transition form factors calculated 
in the nonlocal chiral quark model [50]. I5T] with the mass parameters Mq = 135 MeV and 
strange quark mass M s = 250 well describes the BABAR data. 
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FIG. 6. The pion, r], rf , and ^distribution amplitudes for the nonlocal chiral quark model with parameters 



M u = 135MeV, M s = 250 MeV, M c = 1550MeV, 



From above results one may conclude that the possible origin of the difference of the 
asymptotic behavior of the pion and r] from one side and of the rj' meson form factors is the 
fact that rj' is much heavier than the pion and rj. We have checked that it is this property 
responsible for the change of the almost flat DA in the pion and rj cases for the suppressed 
at end points DA in the cases of the t]',t] c mesons (Fig. [6]) [65] . 
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V. PION TRANS VERSITY FORM FACTORS 



The transversity form factors (TFFs) of the pion provide valuable insight into chirally- 
odd generalized parton distribution functions (GPDs) as well as into the non-trivial spin 
structure of the pion. These interesting quantities have been determined for the first time 
on the lattice [62]. Formally, the TFFs, denoted as B^ ni (t), are defined as 

n-i ^ ^ B K ' u (t'\ 

(7r + (P / )|OT^ 1 '"^ 1 |7r + (P)) =ASP»A V V A m • • • A^P^ +1 • • • P^~ l TnA ' , (10) 

. ' m w 

i=0, even 

where P' and P are the momenta of the pion, the P = \{P' + P), A = P' — P, and t = A 2 . 
The symbol AS denotes symmetrization in u, . . . ,/i n -i, followed by antisymmetrization in 
fi, v, with the additional prescription that traces in all index pairs are subtracted. The 
tensor operators are given by 

where D = ^ denoting the QCD covariant derivative. 

The available full-QCD lattice results [H2] are for B^ u and and for —t reaching 
2.5 GeV 2 , with moderately low values of the pion mass ~ 600 MeV. The calculation 
uses the set of QCDSF/UKQCD Nf = 2 improved Wilson fermion and the Wilson gauge- 
action ensembles. 

Within the effective nonlocal model one gets from the triangle diagram corresponding to 



(10) the transversity pion form factors (see (631 El], where notations are explained) 



An 2 f% i\ (n - I - i)\ J A n + 2 
[2aG m>0fi («, 7) + /3Co, m ,o (a, (3, 7)], (11) 



where i — 0, 2, ... < n — 1 and pion transversity generalized parton distribution (Fig. [7] 

^ (x, e, t) = e (X + — ^ / d 7 / da : e ;f 

0^771,0,0 ( a , ,m,o ( a 5 7) + 7^0,0,771 (a, 

A 2 (l-X) 



(12) 



where £ = , (a 2 = 0), = g±3jg£=fll A = [a + 7 + £ (a — 7)] /(l — X). These 

studies show how the spinless pion acquires a non-trivial spin structure within the framework 
of chiral quark models. 
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VI. CONCLUSIONS 



The main conclusion is that the study of pion and other pseudoscalar mesons provides 
reach information about dynamics of strong interactions and in some rare processes testifies 
the standard model. New high statistics experiments in wider kinematical region are urgent 
for further progress. It is expected new interesting theoretical results of the nonperturbative 
QCD dynamics from the Lattice QCD and effective chiral quark models. 

We are grateful to the Organizers and personally S. Dubnicka, A.-Z. Dubnickova, E. 
Bartos, M. Hnatic, N.Shumeiko, N.Maksimenko, V.Mossolov, S. Eidelman for nice meetings 
and kind invitation to present our results. This work is supported in part by RBFR grants 
10-02-00368 and 11-02-00112. 
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